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COMPARISON OF DECAY OF SOLUTIONS TO TWO 
COMPRESSIBLE APPROXIMATIONS TO NAVIER-STOKES 

EQUATIONS 


CESAR J. NICHE AND MARIA E. SCHONBEK 


Abstract. In this article, we use the decay character of initial data to com¬ 
pare the energy decay rates of solutions to different compressible approxima¬ 
tions to the Navier-Stokes equations. We show that the system having a nonlin¬ 
ear damping term has slower decay than its counterpart with an advection-like 
term. Moreover, me characterize a set of initial data for which the decay of 
the first system is driven by the difference between the full solution and the 
solution to the linear part, while for the second system the linear part provides 
the decay rate. 


1. Introduction 


The Navier-Stokes equations for an incompressible homogeneous fluid in M 3 


d t u + (u-X7)u = Au — X7p, 


div u = 0, 

u o(x) = u(x,0) 


( 1 . 1 ) 


have been extensively studied because of their importance in modelling a wide 
range of phenomena in Fluid Mechanics. Taking divergence in the first line of m, 
using the divergence-free condition and then inverting the Laplacian, we obtain the 
nonlocal relation p = —A~ 1 div (u ■ V)u, which poses very hard problems when 
trying to solve these equations numerically. 

In order to avoid these problems, Temam [II proposed a model approximating 
m in which the pressure and the velocity are related through ep = — divu , for 
e > 0, thus breaking the nonlocality. To “stabilize” this system, i.e. to have an 
energy inequality, he added the nonlinear term ^(div u e )u e , which then leads to the 
compressible system 



( 1 . 2 ) 
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This system has been used in many numerical experiments and has also been the 
subject of some articles concerning its analytical properties (see Fabrie and Galusin- 
ski [2] , Plechac and Sverak 0 ). Recently Rusin [3] proved existence of global weak 
solutions in R 3 and their convergence in L 3 oc (R 3 x R + ), when e goes to zero, to a 
suitable (in the sense of Caffarelli-Kohn-Nirenberg 0 ) solution to the Navier-Stokes 
equation. 

Given the differences in the linear parts and nonlinear terms of (EH) and EH it 
is natural to ask whether these affect the decay rates of the L 2 norm of solutions. 
A useful tool to try to answer this question is the decay character of the initial 
datum, introduced by Bjorland and M.E. Schonbek [6] and refined by Niche and 
M.E. Schonbek |7|. Roughly speaking, the decay character is a number associated to 
every uq £ L 2 (R 3 ) that describes the behaviour of uo near £ = 0, which characterizes 
the norm decay of solutions to linear systems Ut = Cu for a wide class of linear 
operators that includes C = A and C = A + \ V • div. Using the sharp decay rate 
estimates obtained by Niche and M.E. Schonbek and the Fourier Splitting method, 
the following Theorem can be proved. 

Theorem 1.1. (Niche and M.E. Schonbek [TJj Let u be a solution to either it 1.1}) 
or EH), with uq £ L 2 (R 3 ) and decay character r* = r*(uo), with — | < r* < oo. 
Then 


ll u Wlli a (R 3 ) < C(1 + t)- min ff+ r *-§}. 

Remark 1.2. Bjorland and M.E. Schonbek [6] proved this result for the Navier- 
Stokes equations EH- 

The decay rate obtained in Theorem 11.11 provides plenty of information about 
the similarities between eh and E2D- From Theorem 12.101 we have that 

||e t£ uolH 2(R 3) < C(l+f)-(^*), r*=r*(u 0 ), 

so we see that for r* < 1 the linear parts are the ones that have slower decay, 
while for r* > 1 the nonlinear terms are the ones driving the decay due to the 
fast dissipation provided by e tc . Also, the stabilizing term ^(div u e )u e , needed for 
having an energy inequality in EH, does not change the relative strength of the 
linear parts and the nonlinear terms with regards to their influence on the decay of 
energy. We can then conclude that, regarding energy dissipation, the Navier-Stokes 
equations EH and Temam’s approximation EH have the same behaviour. 

The main goal of this article is to study the decay of a different compressible 
approximation to EH and compare its decay rates to those from Theorem 11.11 
using the decay character. More precisely, consider the system 


= Au e -\ — X/-divu e , a > 0 
e 

= Uq(x), 


(1.3) 


d t u e + (it £ • V)u e + a \u € \ 2 u e 
u e (x, 0) 
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introduced by Lelievre 0 and Lemarie-Rieusset and Lelievre :9] as a modifica¬ 
tion of one used by Vishik and Fursikov [TO] to construct statistical solutions to 
Navier-Stokes equations. System ed has the same scaling as Navier-Stokes equa¬ 
tions, thus allowing for its analysis in spaces which contain homogeneous initial 
data leading to selfsimilar solutions. Note that this system differs from Temam’s 
approximation ( 11 . 21 ) in that, instead of having an advection-like term ( divu e )u e , 
it has a strongly nonlinear damping term \u e \ 2 u e . Lelievre [8j proved that for 
Uq £ L 2 (R 3 ), there exists a global in time weak solution to (11.31) and that when a 
and e go to zero, solutions to (11.31) converge, as distributions, to a suitable solution 
to the Navier-Stokes equations. 

We now state our main result which will allow us to compare the behaviour of 
ED and (11.21) to that of (11.31) . 

Theorem 1.3. Let Mq £ L 2 (R 3 ), with r* = r*(uo), with — | < r* < oo. Then for 
any weak solution to (mi we have that for a > C = C(a, e, ||ttolli 2 ) an d an D 
S > 0 


\\u e {t)f L 2 < C(l + t)~ min{ i +r *’i- s} . 

Thus, the strong nonlinear damping |u £ | u e also leads to decay (see the energy 
inequality (13.81) 1. However, there are significant quantitative and qualitative differ¬ 
ences regarding the decay of solutions. First, note that for initial data with r* > 0, 
solutions to Navier-Stokes equations ED and Temam’s approximation ED have 
faster decay rates than solutions to Lemarie-Rieusset and Lelievre’s system. This 
is so precisely due to the presence of this nonlinear damping, which slows down 
decay, as can be explicitly seen in the proof by comparing (13.131) to (13.141) . 

Moreover, for Navier-Stokes equations and Temam’s approximation, the linear 
part determines the decay rates for r* < 1, while for Lemarie-Rieusset and Lelievre’s 
system the linear part is the leading one only for r* < 0. Hence, for Uq with 0 < 
r* < 1 not only the decay rates are different, but so is the dissipation mechanism, 
given by the linear part in the first case, and by the nonlinear terms in the second 
case. In Remark 12.81 we show that Vo = A s Mq, where Uq is in L 1 (R 3 ) D JL S (R 3 ) 
with 0 < s < 1, has decay character r*(v o) = s, thus providing explicit examples of 
initial data that lead to such behaviour. 

Remark 1.4. For any fixed a > 0, when e goes to zero we obtain the Navier-Stokes 
equations with an extra damping term. The decay of solutions to this system has 
been recently addressed, see Cai and Lei E Jia, Zhang and Dong [12], Jiang and 
Zhu [T3] , Jiang [T4]. Theorem 11.31 improves and generalizes the results obtained in 
the case of /? = 3 in [14] . 

This article is organized as follows. In Section[2]we recall some existence results 
and properties of solutions to ED and, following Niche and M.E. Schonbek [7], 
we provide the definitions and results we need concerning the decay character and 
the characterization of decay of linear systems. In Section [3] we prove our main 
Theorem 11.31 
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2. Settings 

2.1. Solutions to (II.3D . For the sake of completeness, we first recall results con¬ 
cerning existence of solutions to (11.31) and their convergence to solutions to (11.11) . 

Theorem 2.1. (Theorem 3.3, Lelievre [8],) Let uo £ L 2 (R 3 ). Then there exists 
a distributional solution u e to i 1 ., *?!) such that u e £ L°°(R +1 L 2 ) DL 2 (R+,Ff 1 ) D 
L 4 (R+,L 4 ). 

Theorem 2.2. (Theorem 4-3, Lelievre [8],) Let uq £ L 2 (R 3 ), with V-ito = 0- Then 
solutions from Theorem 12.11 converge when a and e go to zero, as distributions, to 
solutions to m- 

2.2. Decay character. As the long time behaviour of solutions to many dissipative 
systems is determined by the low frequencies of the solution, Bjorland and M.E. 
Schonbek [6] introduced the idea of decay character of a function uo in L 2 (R”) in 
order to characterize the decay of solutions to Navier-Stokes equations with that 
initial datum. Recently, Niche and M.E. Schonbek [7] generalized this notion in 
order to use data in H s (R n ),s > 0 and to obtain results for other equations. We 
recall now these definitions and results. 

Definition 2.3. Let uq £ L 2 ( R”) and A = (—A)^. For s > 0 the s-decay indicator 
P r s (u 0 ) corresponding to A s u 0 is 

Pf(u 0 ) = limp- 2 ’'-" / |£| 2 s |u 0 (£)| 2 ^ 

Jb(p) 

for r £ (— j + s,oo), where B(p ) is the ball at the origin with radius p. 

Remark 2.4. Setting r = q+s, we see that the s-decay indicator compares | A s uo(£)| 2 
to m = ki 2(9+s) near £ = 0. When s = 0 we recover the definition of decay 
indicator given by Bjorland and M.E. Schonbek [6]. 

Definition 2.5. The decay character of A s uq, denoted by r* — r*(ug ) is the unique 
r £ (— j + s, oo) such that 0 < Pf(uo) < oo, provided that this number exists. 
If such Pf(u 0 ) does not exist, we set r* = — ^ + s, when Pf(uo) = oo for all 
r £ (—f + s, oo) or r* = oo, if Pf(uo) = 0 for all r £ (—f + s, oo). 

Remark 2.6. (Examples 2.5 and 2.6, Niche and M.E. Schonbek |7j). Let uq £ 
L 2 { R") such that u 0 (f) = 0, for |£| < 5, for some S > 0. Then, Pf(uo) = 0, for any 
r £ (— 7- + s, oo) and r*(uo) = oo. If uq £ L P (R”) D L 2 (R"), with 1 < p < 2, then 
r*(uo) = — n ^1 — I'j , so if uq £ L 1 (R") D L 2 (R n ) we have that r*(ug) = 0 and if 
uo £ L 2 (R") but uo $. L P (R"), for any 1 < p < 2, we have that r*(uo) = — §■ 

Let uq £ ff s (R n ), with s > 0. As A s (uo)(£) = |£| s uo(£), the heuristics for the 
decay character given in Remark 1*01 lead us to expect that r*(uo) = r*(A s uo ) = 
s + r*(uo). This is the content of the following Theorem. 

Theorem 2.7. (Theorem 2.11, Niche and M.E. Schonbek [7]J Letuo £ H s ( R”), s > 

0. 
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(1) If —j < r*(u o) < oo then —^ + s < r*(uo) < oo and r*(uo) = s + r*(uo). 

(2) r*(uo ) = oo if and only if r*{uf) = oo. 

(3) r*(uo ) = —\ if and only if r*(uo) = r*(uo) + s = —^ + s. 

Remark 2.8. We can now justify the assertion made after the statement of Theo¬ 
rem 11.31 about initial data for which the decay of solutions to (11.21) and m are 
quantitatively and qualitatively different. Let uq £ L 1 (R. ra ) D L 2 (R") such that 
A s ito £ L 2 (M n ), with 0 < s < 1. From Remark 12.61 we have that r*(uo) = 0, while 
from Theorem rm we have that r*(A s zio) = r*(iio) = s + r*(uo) = s. This proves 
that A s zio has decay character 0 < r* < 1. 

2.3. Linear operators and characterization of decay. We describe now the 
linear operators for which we characterize the decay of the L 2 norm of solutions 
in terms of the decay character. For a Hilbert space X on R", we consider a 
pseudodifferential operator £ : X n —> (^(R"))™, with symbol Ai (£) such that 

(2.4) M(0 = P- 1 (tM0m, Z-a.e. 

where P(f) £ 0(n) and £>(£) = — Ci|£| 2a £,j, for a > c > 0 and 0 < a < 1. Taking 
the Fourier Transform of the linear equation 


(2.5) v t = £v, 

multiplying by v, integrating in space and then using (12.41) we obtain 

^Ht)||| 2 <-c£|^| 2 “|u| 2 ^ 

which is the key inequality for using the Fourier Splitting method in the proofs. 
Remark 2.9. The fractional Laplacian on vector fields in R." 


{£u) i = (-A) a Ui, i = 1, • • ■ n, 

provides an example of (12.41) . as its symbol is (A4(£))j ■ = — C\£\ 2 a Sij. The operator 

Cu = A u H— Xdiv u, e > 0, 

e 

i.e. the linear part of & provides a second example, as (12.41) holds with (A4(^)) 4 • = 
H£l 2 $ij - D{£) = diag(—\£\ 2 ,—\t;\‘ 2 , — (l + \) |£| 2 ) and 


P(f) = 


~&2 

Vf*i 

0 


-Si A 

1 -S 3 


62 

C.3 


where v = (£ 1 ,^ 2 , £ 3 ) has norm one. As a result of this, 
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see Rusin |4|. 

We now state the Theorem that describes decay in terms of the decay character 
for linear operators as in (12.41) . 

Theorem 2.10. (Theorem 2.10, Niche and M.E. Schonhek [Tjj Let vq £ L 2 (R n ) 
have decay character r*(v 0 ) = r*. Let v(t) be a solution to \2.5l) with data Vq- 
Then: 

(1) if < r* < oo, there exist constants Ci, C 2 >0 such that 

C 1 (l+t)-^(5+ r ‘) < ||u(f)||| 2 <C 2 (l + f)-^ +r *); 

(2) if r* = —j, there exists C = C(e) > 0 such that 

\\v{t )\\ 2 L 2 > C(l + t)~ e , Ve > 0, 

i.e. the decay o/||r(t )||| 2 is slower than any uniform algebraic rate; 

(3) if r* = oo, there exists C > 0 such that 

\\v(t )\\ 2 L 2 <C(l+t)~ m , V?7T>0, 

i.e. the decay of Wvft)]]^ is faster than any algebraic rate. 

3. Proof of Theorem 11.31 

Proof: The proof is based on the Fourier Splitting method, introduced by M.E. 
Schonbek to study decay of parabolic conservation laws m and of Navier-Stokes 
equations m , m- As is usual in this context, we prove the estimate assuming 
the solutions are regular enough (for the existence of regular approximations, see 
Lclicvre 0 ). The limiting argument that proves the estimate for weak solutions 
follows that for the Navier-Stokes equations in pages 267-269 in Lemarie-Rieusset 
jl8l and the Appendix in Wiegner m, which we refer to for full details. 

We first show that solutions obey an energy inequality. As 

[ u e {u e -V) u e dx = —— I \u e \ 2 divu e dx, 

Jr 3 2 J R 3 

it follows that 


u e ( u e ■ V) u e dx 


(3.7) 

Multipyling (ESI by u\ 


< i f \u e \ 2 divu e dx 
2 J R3 

< \\\W\ 2 \\LA\ divue \\L 3 < |l|M e |ll4 + ^11 divu e 
integrating in space and using m we obtain 



2 

l 2 - 


(3 - 8) \jt^ {t) ^ - -H Vue WHi 2 - -(«-!) ll« £ Wlli*. 

which, given the hypotheses on e and a, provides an energy inequality. Then 


~\\^mi 2 <-\\^mi 2 


C\\div u e \\ 2 L 2 ■ 


(3.9) 






DECAY OF SOLUTIONS FOR APPROXIMATIONS TO THE NSE 


7 


Now let 


Bit) = {{ 6 K 3 : l£P < jgly} 

where r is a positive increasing function with r(0) = 1 and C is an appropiate 
constant. From (13.91) we have that 

li\\u e {t)\\ 2 L 2 <-C [ [ \u*(U)\ 2 dti. 

1 dt J R3 Zr \t) J B(t) c 

Adding and substracting a term similar to the one on the right side of this inequality, 
only that with B(t) as the domain of integration, and then mutiplying by r(t) we 
obtain 

(3.10) ^Ar(t)\\u*(t)\\ L *)<r'(t) [ |£(£,i)| 2 d£. 

at J B(t) 

We now prove a pointwise estimate for 

u e (£,t) = e tM ^u\ q(^) — f s ) ds 

Jo 

where e tM ^ is as in (EH) and 

G(£,s) = V)u e + |« £ |V) , 

where J- is also the Fourier transform. As 


( u e ■ X7)u e = V ■ (u e <8 u e ) — ( divu e )u € 

and 


we obtain 


(div u e )u e (t;, t) 


< \\divu e {t)\\ L 2 \\u e (t)\\ L 2 , 


■ v)« e (t)) i < + l -\\dwu\t)\\ L 4u\t)\\ L , < 

Also 


T 


(V|V) (t )| < ||-F(V|V) (t)|U~ < II K(t)|V(t)|| L i < IKWIlla. 


Now, we estimate the nonlinear transport term. We have 


f e {t ~ s)M ^T (( u e ■ V)u e )(^s)ds < cf 

Jo Jo 


o-C(t-8) ICI 2 


l£IK(s)lli» ds 


(3.11) 


< c\z\(j*\\u*(s)\\i 2 dsy 


Suppose now that 
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(3.12) \\u e {t)\\ 2 L2 < C(l + t)~ p , 

for some /3 > 0 and /3 ^ 1. We then have, after choosing r(t) = (t + l) 3 

/B( t) (f e(MA1K)j((,i£ ' vK)(e ’ s)d ') d( - c i^M ) 2 {i+t)2{l ~ 0) 

(3.13) < C{1 

where we used (13.111) and (13.121) . We now estimate the damping term. We have 


X( ) (/ g(t (£.s)ds) d£<C J (^J o KWIIi-*) d £ 



< C j Qf* \\u^s)\\ 2 l2 ds^j (j* ||u £ (s)||l4 ds^j dZ 

(3.14) < C(t + 1)“(^ +/3 ) 

where we used the interpolation 

||u £ (t)|U3<|| M £ (t)||l 2 || U £ (t)||| 4 , 

Holder’s inequality, (13.121) and the fact that u e £ L^°L 2 nLfL^.. Of these two terms, 
(13.141) has the slower decay. 

The only apriori estimate we have is ||u e (t )||^2 < C, i.e. /3 = 0. So, in the ball 
B(t) we have 



/ 

|W £ (£,*)| 2 ^ < 

cf \e tM ^\ 2 dC 


J B(t) 


J B(t ) 



+ 

C [ ( [ e (t - a)M( ‘ } G(€,s)ds) 




JB(t) \J 0 / 

(3.15) 


< 

C , (i + l)-(* +r ’) + C(t + l)-3, 

which leads to 



(3.16) 


\\u*m 2 L> 

< C(1 +t)“ min{ i +r *^ } . 


di 


We now use this to boostrap our decay estimates by improving the value of /?. 
When r* < —1, we have that /3 = § + r* and (13.141) provides no improvement. 
When r* > —1, then B = \ which leads to 
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(3.17) 

Hence 



< C f \e tM{i) ud\ 2 d£ 

JB(t) 

+ c [ ([ e( t - s)M ^ ) G{Z,s)ds) 

J B(t) \J 0 J 

< C(t + l)-^ +r ') + C(t + l)~ 1 . 


(3.18) \\u e (t)\\ 2 L2 < C(l + t)- min {|+ r *’ 1 >. 

Again, r* < — i does not produce an improvement on the decay while r* > — i 
leads to (3 = 1, for which (13.1311 and (13.141) do not hold, as 


/ ||w e (s)|li 2 ds < Cln(t + 1). 

Jo 

However, working along the lines of these estimates we obtain 


(3.19) / B ( t )(( eMMK)j (( ue ' V )“ £) (^ )(is ) d ^- C (^B) 2 ln 2( * + l) 

and 

t 2 3 . 

(3.20) J ( jf ( \uf (£, s) ds ) d£ < C * ln(i + 1). 

Again, the second term has slower decay and choosing r(t) = (t + l) 7 for large 
enough 7 > 0 we have that, as ln(t + !)<(£ + I) 5 for any S > 0 



/ \u^,t)\ 2 d£ < 

C / |e 

■ tM ^\ 2 dt 


J B(t) 

J B(t) 



+ 

c f ( 

f e< t -') M(i) G(€,s)ds) d£ 



J B(t) \ 

Jo J 

(3.21) 

< 

c(t + iy 

-(i+ r *) + c(t + lyd- 6 ). 

Then 





\\u £ (t)\\ 2 L2 <C(l+t)- m ' ln{ % +r *’i~ s \ Vd > 0. 

This boostrapping process does not improve the decay further, as now for r* > 0, 
the integral of ||u £ (t)||i 2 is bounded by a constant. □ 
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